Laser cooling of a trapped particle with increased Rabi frequencies 
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This paper analyses the cooling of a single particle in a harmonic trap with red-detuned laser 
light with fewer approximations than previously done in the literature. We avoid the adiabatic 
elimination of the excited atomic state but are still interested in Lamb-Dicke parameters n <C 1. 
Our results show that the Rabi frequency of the cooling laser can be chosen higher than previously 
assumed, thereby increasing the effective cooling rate but not affecting the final outcome of the 
cooling process. Since laser cooling is already a well established experimental technique, the main 
aim of this paper is to present a model which can be extended to more complex scenarios, like 
cavity-mediated laser cooling. 

PACS numbers: 03.67.-a, 42.50.Lc 



I. INTRODUCTION 

The idea of laser cooling for neutral atoms was first 
suggested by Hansen and Schalow [l| and independently 
for trapped ions by Wineland and Dehmelt 0]. Since 
light carries momentum as well as energy, scattering light 
on matter can result in significant changes of the vibra- 
tional energy of massive particles. Over the last decades, 
laser cooling techniques have been designed which al- 
low the cooling of atoms and ions to the the micro and 
nanokelvin temperatures needed for quantum coherence 
and degeneracy Examples of which are Sisyphus 

cooling [(| and evaporative cooling Q. In this paper we 
focus our attention on laser sideband cooling of a particle 
in a harmonic potential Q which can be used for exam- 
ple to transfer single ions to their motional ground state 
with a very high precision [t| ■ 

The idea of laser sideband cooling of a single trapped 
particle is relatively straightforward. If the trapping po- 
tential of the particle is strong enough for its motion to 
become quantised and if the corresponding phonon fre- 
quency is much larger than the spontaneous decay rate 
of the excited electronic state, it becomes possible to re- 
solve the motional sidebands with the cooling laser. As a 
consequence, the excitation of the excited electronic state 
of the trapped particle on the red sideband is most likely 
accompanied by the annihilation of a phonon. When 
followed by the spontaneous emission of a photon, the 
particle returns into its electronic ground state without 
regaining a phonon which implies cooling. The cooling 
cycle only stops when the particle reaches a state with 
almost no phonons. The final population of excited vi- 
brational states is in general very small and only due to 
highly off-resonant excitations of the ground state of the 
atom-phonon system. 

The theory of laser cooling has already been studied 
in great detail in the literature. For reviews on this topic 
see for example Refs. [TTl - flij . The first theoretical dis- 
cussion of laser cooling with red-detuned light based on a 
combination of simple classical and quantum ideas can be 
found in Ref. Q by Wineland and Itano. Lindberg and 
Stenholm later introduced the tool for a full quantum 



treatment of laser cooling by deriving a master equation 
for spontaneously emitting atoms with atomic recoil in- 
cluded [HI (cf. also Refs. (l6l - [20j ). An alternative but 
consistent analysis of the laser cooling of trapped ions in 
a running and in a standing wave configuration has been 
presented by Cirac et al. in Ref. [2l|. The main result 
of these papers is a cooling equation of the form 



- A + ) m + A + 



(1) 



where m denotes the mean phonon number. The A± can 
be interpreted as transition rates between states with dif- 
ferent phonon numbers and relate to actual cooling and 
heating rates. Moreover, one can show that the station- 
ary state phonon rate m ss is given by [HI, BIS III 
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Experimental results confirm the general dependence of 
this stationary state phonon number on the emission rate 
of the excited electronic state of the trapped particle T 
and on its phonon frequency v [23| . 

The purpose of this paper is to analyse the cooling 
process of a single trapped particle in a harmonic po- 
tential with fewer approximations than previously done 
in the literature. For example, we avoid the adiabatic 
elimination of the excited atomic state, thereby obtain- 
ing an analysis of the cooling process which applies for 
a much wider range of Rabi frequencies of the cool- 
ing laser. As a result we find that does not have to 
be small compared to the spontaneous decay rate T and 
the laser detuning A. For example, in the weak confine- 
ment regime with y < T, the Rabi frequency 51 can be 
as large as 0.3 T and in the strong confinement regime 
with r < the Rabi frequency can be as large as 
0.3 v without noticeably increasing the stationary state 
phonon number to ss . This is an interesting observation, 
since the effective cooling rate j c scales as f2 2 which in- 
creases rapidly when Q increases. 

Our calculations are more straightforward than pre- 
vious calculations, since we replace the atomic raising 
operator <j + and the phonon annihilation operator b by 
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two new operators x and y. These correspond to par- 
ticles that are neither atoms nor phonons and commute 
with each other. Most importantly, they provide a rep- 
resentation of the Hamiltonian which no longer contains 
atomic displacement operators. Instead it depends on 



which take non-linear effects in 



terms like x^x(y 
the atom-phonon interaction into account [22j . In the 
following we use this Hamiltonian to obtain a manage- 
able set of cooling equations which are differential equa- 
tions for the time derivatives of expectation values. As in 
Refs. [HI [3 [IE OS [III , we are interested in the dynam- 
ics of the cooling process on the very slow time scale given 
by the cooling rate j c which scales as rj 2 with rj <§C 1 . The 
only approximation involved in the following calculations 
of stationary state phonon numbers and effective cooling 
rates is to neglect higher order terms in the Lamb-Dickc 
parameter rj. 

The main purpose of this paper is to establish and test 
a framework for the modelling of laser cooling which can 
be extended relatively easily to more com plex cooling 
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scenarios like cavity-mediated laser cooling 
study of possible quantum optical heating mechanisms 
in sonoluminescence experiments [25j |. Cavity- mediated 
laser cooling for example has many similarities with laser 
cooling but due to being more complex there is more flex- 
ibility in choosing different scenarios and experimental 
parameters [2a |. 

There are six sections in this paper. Section |TT] intro- 
duces the master equation of a single laser-driven trapped 
particle. Section IIIII uses this master equation to derive 
a closed set of 23 cooling equations. These simplify to a 
set of five equations in the weak confinement regime and 
to a single equation in the strong confinement regime, 
respectively. In Section ITVl we show that the phonon co- 
herences and the mean phonon number m always reach 
their stationary state. Section [V] presents cooling rates 
and stationary state phonon numbers and compares ana- 
lytical and numerical results. Finally, we summarise our 
findings in Section PVTl Mathematical details are confined 
to Apps.|A][D] 



II. THEORETICAL MODEL 

Let us start by introducing the theoretical model and 
the experimental setup. It consists of a single cooling 
laser which drives a strongly confined single particle. As 
long as the trapping potential is approximately harmonic, 
we can describe the motional states of the atom by a har- 
monic oscillator Hamiltonian. In this section, we consider 
the motion of the trapped particle as quantised and in- 
troduce the model which allows us to predict the time 
evolution of its mean phonon number m. The purpose of 
the cooling laser is to minimise the number of phonons 
in the motion of the particle in the direction of the laser. 
Cooling other vibrational modes requires additional cool- 
ing lasers. 



A. The Hamiltonian 

The Hamiltonian of a single trapped particle inside the 
free radiation field and with external laser driving can be 
written as 

H = i?nuclei + ^electron + Afield + #dip ■ (3) 

The first two terms are the free energy of the electronic 
states and of the quantised vibrational modes of the 
trapped particle. The third term describes the energy 
of the surrounding free radiation field. The last terms 
take the dipole interaction of the electronic states with 
the present electromagnetic fields, i.e. the laser and the 
free radiation field, into account. We now have a closer 
look at every term in this equation. 

Suppose the particle is effectively a two-level system 
with ground state |0) and excited state |1) and the en- 
ergies Hloq, hu 1 and Hujk denote the energy of a single 
atomic excitation, of a single phonon excitation, and of a 
single excitation of the modes of the free radiation field, 
respectively. Then we have 



H, 



electron 



H, 



nuclei 
fffield 



hujQ a + a 



(4) 



kA 



where the operators a~ = |0)(1| and a + = |1)(0| are 
the atomic lowering and raising operator and where b 
and dkA are the phonon annihilation operator and the 
annihilation operator of a photon with wavevector k and 
polarisation A, respectively. These operators obey the 
commutator relation 

[«kA,4 A ] = [6,6+] = I (5) 

which is the usual commutator relation for bosonic an- 
nihilation operators. All other photon commutators are 
equal to zero. 

The final term in Eq. ^ describes the dipole inter- 
action of the electronic states |0) and 1 1} of the trapped 
particle with the free radiation field and the applied laser 
field. Within the usual dipole approximation [30(, it can 
be written as 

H dip (t) = eD [E fic i d (R) + E L (R,t)] . (6) 

Here e is the electron charge, D is the dipole moment of 
the particle, i.e. the position operator of its outer electron 
with respect to the atomic nuclei at position R, while 
E n eid(R) and El(R, t) denote the electric field amplitude 
of the free radiation field and of the laser field at time t, 
respectively. The dipole moment D equals 



D = Dm a + H.c. , 



(7) 



where Doi is a 3-dimensional complex vector. The elec- 
tric field operators are given by 



kA 



E fi cid(R) = i^/7^^ekAa k Ae ik - R + 



2e L 3 



El(R, t) = E e i(kL - R - WL *' +c.c. 



(8) 
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with L 3 being the quantisation volume of the free radia- 
tion field and e^x being a unit length polarisation vector 
orthogonal to k. Moreover, E , kL, and wl are the am- 
plitude, the wave vector of length k^, and the frequency 
of the applied laser field, respectively. 



This equation shows that the cooling laser establishes a 
coupling between the electronic states |0) and 1 1) of the 
trapped particle and its quantised motion. As we shall 
see below, the coupling to the free radiation field is the 
origin of spontaneous emission and recoil effects. 



B. The displacement operator 



C. Interaction picture 



In the following, we assume that the incoming laser 
field has the same direction as the quantised motion of 
the trapped particle, since this maximises the effect of the 
cooling laser. Considering this motion as quantised with 
the phonon annihilation operator b from above, hence 
implies 



k T R 



r, (b + 



(9) 



where the Lamb-Dicke parameter rj is a measure for the 
steepness of the effective trapping potential seen by the 
ion [11|. Notice that Eq. © applies as long as the trap- 
ping potential seen by the atom does not depend on its 
respective electronic state. This means, the following cal- 
culations apply to a trapped ion and to a trapped neutral 
atom with a magical wavelength . 

Taking Eq. © into account, we find that the laser 
Hamiltonian is a function of the particle displacement 
operator [28| 



D(ir)) = e - i " (fo+bt) 
which is a unitary operator with 

D(ir])bD(iri^ = b 



D(ir])HD(iri) = b- 



"7, 

IT}. 



(10) 



(11) 



Using the polar coordinates i? and p, putting the z-axis in 
the direction of the cooling laser, and writing the general 
wave vector k of length k as 



we find that 

k R 



sin $ cos <p 
k = k I sin -d sin ip 

COST? 



k sin $ [R x cos tp + R y sin ip] 
rjk cosd 



(12) 



kh 



(6 + 6*) 



(13) 



where R x and R y are the x and the y component of the 
vector R. Writing the Hamiltonian TJdip hi Eq. as a 
function of displacement operators, it becomes 



ff d i p (t) = e [Doitf- +H.c] 



kA 



E* D(ir,) e iWLt 
vqk cos d 



2e L3 



xe 



ik sin $[R X cos tp-\-R y sin <p] 



H.c. (14) 



Before continuing our derivation of the master equa- 
tion, it is convenient to transform the Hamiltonian H in 
Eq. ([3]) into an interaction picture. To do so, we choose 

H a = Huj L <7 + <J~ + iffield (15) 

with -f/fieid as in Eq. (@|. Neglecting relatively fast oscil- 
lating terms as part of the usual rotating wave approxi- 
mation, the interaction Hamiltonian Hi, 



H! = U^t,0) (H-H )Uo(t,0), 



(16) 



becomes 



E t /'i^fccosi? 

kA 

v, „— ife sin i?[ii x cos tp+Ry sin 1/3] Uuik— UL>t 

Ac C 



1 



Hfl D{ir])a- + H.c. + HAa+a^ + hv b f b . (17) 



Here A denotes the detuning between the laser and the 
relevant atomic transition and £1 and g^x, 

2e D i ■ Eq 



.9kA = -ie 



2he L 3 



Doi • £ kA 



(18) 



are the usual laser Rabi frequency and the atom-field 
coupling constant. 



D. Spontaneous emission and recoil 

The interaction Hamiltonian Hi in Eq. (|17p is now the 
starting point for the usual derivation of the master equa- 
tions. Suppose the state of the laser-driven trapped par- 
ticle is at t = given by the density matrix p, while the 
free radiation field is in its vacuum state |0). Taking this 
into account, the density matrix p(Ai) of the particle at 
time At can be written as [2!| 

p(At) = U cond (At,0)pUl nd (At 1 0) + n(p)At, (19) 

where At denotes the typical response time of the envi- 
ronment, i.e. the typical time it takes the environment to 
absorb a photon from the free radiation field, and where 

U cond (At,0) = (0|C/i(A*,0)|0), 

U{p) = il I 5oAl^ <lkA|(7l(A< ' 0)p 



kA 



^|0)(0|i7 I t (At,0)|l kA ). (20) 
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The first term in Eq. (jT9|) describes the subensemble with 
no photon emission in At. The second term in this equa- 
tion is the unnormalised state of the subensemble with 
an emission in At [2jj. Taking the time derivation of 
p(At) on the coarse grained time scale At into account, 
we obtain the usual master equation 



H condP - P Hl nd +K(p), (21) 



in Lindblad form, where H con d is a non-Hermitian con- 
ditional Hamiltonian with U con d being the corresponding 
no-photon time evolution operator. 

In order to calculate H con d and ~R-(p), one usually ex- 
ploits second order perturbation theory. Since the dis- 
placement operator D(irj) is a unitary operator, i.e. 



L>0/)D(if7) t = Diirj)^ D(irj) 



1 



(22) 



the derivation of the conditional Hamiltonian H con( \ re- 
mains exactly the same as in the case, where the motion 
of the particle is not quantised. This means, we find that 



H cond = -W,D(iT])a + H.c. + ftAer+cr~ 



(23) 



where the spontaneous decay rate T of the excited elec- 
tronic state |1) is given by 



_ e u 2 
r - 3^? |Do11 



(24) 



However, the reset operator IZ(p) now contains recoil 
terms. Proceeding as described in App.0 we find that 

Op pi 

K(p) = — J dta-D(ir ) ()pD(ir ] O i <T + 

x [l + |d 3 | 2 + (l-3|d 3 | 2 )C 2 ] , (25) 

where d% denotes the z-componcnt of the normalised 
dipole vector Doi/|Doi|. The above reset operator 
is different from the one often used in the literature 
pH li~9l . l2~fl ] . The reason for this is that the authors of 
these references only consider the case where GJ3 =0, as 
it applies to certain atomic level schemes and laser con- 
figurations. The above reset operator IZ(p) is more gen- 
eral. It is also consistent with the reset operator of a free 
particle whose motion is not quantised. In this case, the 
displacement operator D(irj cos 1?) becomes a number and 
the integration over £ results indeed in 1Z(p) =Y o~ pa + . 
This expression is independent of as it should. 

III. COOLING EQUATIONS 

The master equations which we derived in the previous 
section can now be used to derive differential equations 



for expectation values, so-called rate or cooling equa- 
tions. However, this is not a straightforward task due to 
the presence of the displacement operator D in Eq. (|2"3"|) . 
To overcome this problem, we now introduce two new 
operators x and y which replace the particle and the 
phonon operators <r~ and 6, respectively. Both opera- 
tors describe neither electronic excitations nor phonons 
but provide nevertheless a natural description of trapped 
particles. 



Transformation of the Hamiltonian 



To simplify the Hamiltonian H\ in Eq. (|14j) . we now 
introduce a new annihilation operator x as 



x = D{yq) a 



(26) 



Using the commutator relation ([5]) and the unitarity of 
D in Eq. (gTJ), i.e. the fact that 

D(i7?)D(i77)t = D(i v )D(i v y = 1 , (27) 

one can show that x obeys the commutator relation 

= 1 - 2x t cc. (28) 

The operator x differs from a~ by the fact that its appli- 
cation not only transforms |1) into |0) but also induces a 
kick which displaces the particle. 

Eqs. <j5j) , and (|28|) can now be used to derive for 
example the commutator relations 

[x, b] = — [x, 6 T ] = iri x , 
[ x \b] = - [x t ,6 t ] = -vqx* . (29) 

These can then be used to show that 

[x, b^b] = —b]x(b — b^) — rj 2 x , 

[x\ tfb] = h](b - 6 t )x t + if x^ , 

[x*x,b] = [x*x,tf] = [x*x,tfb]=0. (30) 

The operators x and b and functions of them do not 
commute in general. This means, although substitut- 
ing Eq. (|28|l into the interaction Hamiltonian Hi into 
Eq. (fT4")) , simplifies it a lot, we do not yet have a Hamil- 
tonian which can be analysed easily. 

To overcome this problem we now introduce another 
operator y as 



b — irj x^x . 



(31) 



Using the commutator relations in Eq. (|30[) , one can show 
that y is a bosonic operator which obeys the commutator 
relation 



[y,y f ] = 1- 

Moreover, we now have 

[x,y] = [x\y] = 



(32) 



(33) 
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which can be checked using the commutator relations in 
Eqs. (J2SJ and ([55]). 

Using the notation introduced in this section, the con- 
ditional Hamiltonian H conc i in Eq. (|23|) and the reset op- 
erator TZ(p) in Eq. (f25j) become 

#cond = + xt ) - 'Arjvx { x{y - y 1 ) 

+h (A + ifi/^J a^cc + /i^ y^y — — HT x^x , 

op pi 

TZ(p) = — y dCxD(i»j(l- 0)^-0(^(1-0)^ 
x [l + |rf 3 | 2 +(l-3M 3 | 2 )C 2 ] ■ (34) 

The new conditional Hamiltonian -ff C ond no longer con- 
tains any exponential terms. 

B. Time evolution of expectation values 

The time derivative of the expectation value of a time- 
independent operator A equals 

(A) = Tt(Ap). (35) 

The master equation in Eq. ([21]) hence implies that 

(A) = -~(AH cond -Hl nd A) 

op pi 

+ Y J i d C (^D(i v (l - 0) A D(i v (l - ()?x) 

x [l + |d 3 | 2 + (l-3M 3 | 2 )C 2 ] (36) 

with -f/cond as in Eq. (|34[) . This equation describes the 
time evolution of the expectation value {A) within the in- 
teraction picture which we introduced previously in Sec- 
tion IE] 

In the following we are especially interested in the time 
evolution of the mean phonon number m, 

m = (6+6) . (37) 

Using Eqs. (f2"5f and (pTTj). we find that this expression is 
the same as 

m = n 2 - r] k 12 + if ni (38) 

with ni, ri2, and fci2 defined as 

ni = (x^x), n 2 = (y^y) , fc i2 = i (a: t a;(j/ - ?/)) . (39) 

To predict the time evolution of m, we need to evaluate 
the time evolution of these three expectation values. As 
we shall see below, we only obtain a closed set of cool- 
ing equations, if we consider in addition the expectation 
values 

k 7 = (y + y^), k 8 =i(y-ij), k 9 = (if + j/ t2 ) , 
fc 10 = i(y 2 -y t2 >, k n = {x^x(y + y^} (40) 



and the expectation values defined in App. [Bj Since all 
of these variables are expectation values of Hermitian 
operators, they are real and their time evolution is given 
by real differential equations. 

Let us first have a look at the y operator expectation 
values 77,2 and kj to fcio. Using the master equation in 
Eq. (|5o]l to calculate the time derivatives of these five 
variables, we find that 

h 2 = rjv fcn - ?]T k u + r] 6T ni , 

kj = 2r\v n\ — v kg , 

kg = v ky — 2rjT n\ , 

kg = — 2v fcio + 2r\v fcn + 2i]T k\ 2 — 2rj 2 8Tn\ , 

fcio = 2v k 9 + 2r\v k 12 - 2r]T k n . (41) 

The factor 8 in this equation, 

= "Mai 2 ), (42) 

depends explicitly on the direction of the emitting dipole 
moment. It relates to the parameter a used in previous 
papers [IJ (3 Esi H3 v i & the equation 

9 = l + a~l\d 3 \ 2 . (43) 
5 

The time derivatives of the relevant x and mixed operator 
expectation values can be found in App. [Cj Below we 
simplify these equations with the help of an adiabatic 
elimination of all relatively fast evolving variables. 

C. Weak confinement regime 

We now have a closer look at the case, where the 
trapped particle experiences a relatively weak trapping 
potential and where the Lamb-Dicke parameter r\ is much 
smaller than one. More concretely we assume in the fol- 
lowing that 

K<T and 77 < 1 . (44) 

In addition we assume that the Rabi frequency f2 and 
the detuning A are at most comparable to T and defi- 
nitely not much larger. However, notice that we do not 
demand that f2 is much smaller than T. A closer look at 
the cooling equations in App. [Cl shows that this choice of 
parameters causes the y operator expectation values 77-2 
and kj to fcio to evolve on a much slower time scale than 
all other relevant expectation values. The reason for this 
is that these variables are all x or mixed operator expec- 
tation values which decay with the spontaneous atomic 
decay rate V. 

Taking this into account and eliminating 774, k\, k 2l 
and fci 3 to /C24 adiabatically from the system dynamics, 
we obtain a closed set of five effective cooling equations 
which applies after a relatively short transition time and 
which can be written as 

{n 2 ,k 7 ,k 8 ,k 9 ,k 10 ) = M (n 2 , fc 7 , k 8} kg, k w ) 

+ (/? 1 ,/3 2 ,/?3,/34,/3 5 ) T • (45) 
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A closer look at Eq. (|4Tj) shows that the time derivatives 
of the y operator expectation values n 2 , and k 7 to fcio 
depend only on m, fcn, and fci2- The calculation of the 
5x5 matrix M therefore only requires the calculation of 
rii, fen, and fci2 which can be found in App. [C] Substi- 
tuting Eqs. Ip5)l . JDS]), (ICl|, and ([UTTjl into Eq. (|4T|) . 
we find that A/ can be written as 



M 



/ 


(2) 
«11 


"12 


"13 


(2) 





\ 










— V 















V 


(2) 
"33 












(2) 
"41 


"42 


"43 


(2) 






V 





a (1) 

"52 


a (1) 

"53 


2v 


(2) 
"55 


/ 



(46) 



The first order matrix elements a\- in this equations are 
given by 



(i) (i) (i) (i) 

"52 = -«43 . "53 = "42 



(47) 



with /i 2 defined as in Eq. (|C3[) . The non-zero matrix 

(2) 

elements a;, of M in second order in 77 and in first order 
in ^ are given by 



(2) _ (2) _ (2) _ (2) _ 1677 2 i/Ar» 2 
"ll — "33 — "44 ~ "55 



2) 8?7 2 zyArfi 2 (2) 32r; 2 ^ArO 



(2) 

l 14 — "4 ! "41 

M 



(48) 



D. Strong confinement regime 

Let us now have a closer look at the parameter regime 
where the phonon frequency v and the detuning A exceed 
the spontaneous decay rate T and the Rabi frequency 
by at least one order of magnitude, 



fi, T < j/, A , while 77 < 1 



(51) 



In this so-called strong confinement regime, the time 
scale separation in the dynamics of the trapped particle 
is different than in the previous subsection. The adia- 
batic elimination performed in the previous section does 
not apply. However, at least at the end of the cooling 
process when 712 is already very small, we can assume 
that the expectation values rii, TI4, fci, k 2 , and k? to k 2 4 
evolve much faster than the y operator population n 2 . 
This means, we can simplify the system dynamics via an 
adiabatic elimination of all expectation values other than 
7i2- Doing so (cf. App.|D]), we obtain the effective cooling 
equation 

n-2 = -7c n 2 + c. (52) 
The frequencies "f c and c in this equation are given by 

v 2 m 2 rfm 2 



7c 



4(A-i/) 2 4(A + ^) 2 
77 2 m 2 



4A 2 



A 2 



(A + v) 



1 



(53) 



up to second order in r\. It is obvious that 7 C is the 
effective cooling rate for strongly confined particles. 



Proceeding as above, one can show that fi\ — p\ is up 
to second order in 77 to a very good approximation given 
by 

M 

Moreover, we find that the coefficients f3 2 to /3s in Eq. (|45|) 
equal 

(i) _ 2 W n 2 ( i)_ 277m 2 



ft™ = pP = (50) 

in first order in 77. We now have a closed set of five differ- 
ential equations which can be used to analyse the time 
evolution of the y operator expectation values analyti- 
cally and numerically. Notice that the weak confinement 
regime which we introduced in Eq. ((44]) does not allow 
for the adiabatic elimination of the y operator coherences 
k 7 to fcio, since these evolve in general on the same time 
scale as the y particle population n 2 . 



IV. STABILITY ANALYSIS 



In the strong confinement regime, the cooling process 
is governed by a single differential equation (cf. Eq. ([5"2"j) ) 
with an always positive cooling rate 7 C and an always pos- 
itive stationary state phonon number m ss . This means, 
in the strong confinement regime, the trapped particle 
always reaches its stationary state. However, it is not 
clear whether the same applies in the weak confinement 
regime, where the cooling process is described by five lin- 
ear differential equations (cf. Eq. (|4"5)) ). In this section, 
we therefore have a closer look at the dynamics induced 
by these equations. 

To do so, we introduce the shifted y operator expecta- 
tion values 

(n 2 , fc 7 , fc 8 , fc 9 , fci ) T = (n 2 ,k 7 ,k$,k 9 ,k 10 ) T 

+M- 1 (/3i,/3 2 ,/3 3 ,/34,&) T . (54) 

This definition means that the tilde and the non-tilde 
variables differ only by the stationary state solutions of 



7 




FIG. 1: Diagrams illustrating the time evolution of the expectation values kj to fcio, and hi for A = 0.5 T, v — 0.1 V, 
— 0.01 F, and ds = 0. It is assumed that the y particles are initially in a coherent state. In (a), only terms in zeroth 
order in 77 have been taken into account. As expected, we find that the mean phonon number remains constant in time. 
The coherences kt to fcio evolve such that their points in the respective phase diagrams lie on circles. This means that 
the y particles remain in a coherent state. In (b), only terms in first order in 77 have been taken into account. All five 
eigenvalues of M are still either zero or purely imaginary which is why there is still no reduction of 7I2. In (c), also the 
second order terms in 77 are taken into account. The coherences hy to kio now evolve towards zero. We now observe an 
exponential decrease of 712- This implies a reduction of the mean phonon number m in zeroth order in 77, ie. cooling. 



the non-tilde expectation values. Substituting Eq. 
into Eq. 1(15]). we find that 



(h 2 ,kr,k 8 ,kg,k 



10 J 



M 



(™2, 



kr, k 8 , k 9 , ki 



(55) 



The stationary state solution of these effective cooling 
equations is the trivial one with all variables equal to 
zero. However, notice that the corresponding stationary 
state is only reached, if all eigenvalues of M have negative 
eigenvalues. 



A. Time evolution for 77 = 

We first have a closer look at the time evolution of the 
y operator expectation values for 77 = 0. One can easily 
check that the eigenvalues of the matrix M in Eq. (|46p 



are in this case simply given by 

Ai = 0, A 2 , 3 = Ti^, A 45 = qF2if. 



(56) 



Taking this into account and solving Eq. (|55[) analyti- 
cally, we find that 



na(*) 
~k 7 (t) 

h(t) 

kg(t) 



MO), 

cos vt — sin vt 
sin vt cos vt 



cos 2vt 
sin 2vt 



- sin 2vt 
cos 2vt 



k 7 (0) 

*b(o) 
*b(o) 

fcio(0) 



(57) 



This behaviour is illustrated in Fig. HJa) which shows 
numerical solutions of the effective cooling equations in 
Eq. (|55[) for the case where the y particles are initially 
in a coherent state. The first two phase diagrams show 
ks and kio as functions of kj and kg, respectively. The 
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fact that all points lie on a circle means that the state of 
the y particles remains (at least approximately) coherent 
throughout the cooling process. Fig.QJa) moreover shows 
that the y particle population hi and therefore also the 
mean phonon number m (cf. Eq. ([38[) ) in zeroth order in 
r\ remains constant in time. This is exactly as one would 
expect. There cannot be any cooling without a coupling 
between the electronic and the vibrational states of the 
trapped particle. Higher order corrections in r\ have to 
be taken into account. 



B. First order corrections 

Calculating the eigenvalues of M in Eq. (|46|) up to first 
order in 77, we obtain again Eq. ([55)). All of them are ei- 
ther zero or imaginary. However, there are some small 
corrections to the eigenvectors in first order in 77. As a 
result, the shifted y particle population hi remains no 
longer constant in time but oscillates on the time scale 
given by the phonon frequency v. This is illustrated in 
Fig-QTb) which shows a numerical solution of the effective 
cooling equations in Eq. (|55j) with all first order correc- 
tions in 77 taken into account. However, since the eigen- 
values of M have no real parts, hi does not reach its 
stationary state value. No cooling occurs. 



C. Second order corrections 

Calculating again the eigenvalues of the matrix M in 
Eq. ([46]) but now taking also the terms in second order 
in r\ into account, we find that 

\ (2) 

Ai = a\{ , 
A 2 ,a = ^T^-(af>) 2 , 

A 4 , 5 = af 1 ) TiV // 4^ 2 -ai4 ) aii ) - (58) 
(2) 

Since the matrix element <x\{ i s always negative, all 

five eigenvalues of M have negative real parts. This 

means, all tilde variables are damped away and tend to 

(2) 

zero on the time scale given by l/a\i ■ This is illus- 
trated in Fig. []Jc) which illustrates a numerical solution 
of Eq. ([55]). 

Since the y coherences kj to fcio do not increase in time 
but oscillate instead with a slowly decreasing amplitude 
around constant values, the cooling process remains sta- 
ble and the trapped particle eventually reaches its sta- 
tionary state. This observation is taken into account in 
the following section, where we analyse the cooling pro- 
cess in more detail. It allows us to replace the coherences 
kj to fcio by their time averages. This means, in the weak 
confinement regime, the calculations in the following sec- 
tion apply only towards the end of the cooling process, 

(2) 

i.e. after a transition time of the order ofl/a^ . 



V. COOLING RATES AND PHONON 
NUMBERS 



Taking the results of the previous section into account, 
replacing the y operator coherences k*t to fcio by their 
time averages, and adiabatically eliminating all rapidly 
evolving expectation values, we obtain an effective cool- 
ing equation of the same form as the effective cooling 
equation in Eq. ([52]) . The reason for this is that the time 
averages of the coherences ki to fcio are exactly the same 
as the quasi-stationary state solutions obtained when set- 
ting their time derivatives in Eq. (|4"5j) equal to zero. The 
solution of Eq. (|52p implies that the mean phonon num- 
ber at time t can to a very good approximation be written 
as 

m(t) = [m(0) - to ss ] e" 7c * + to ss (59) 

in zeroth order in r\ with to ss = c/^ c , since ni and m 

are the same for 77 1. Notice that the general solution 

in Eq. ([59]) applies in the strong as well as in the weak 

confinement regime, although in the weak confinement 

( 2) 

regime only after a transition time of the order of \/a\-{ . 
In this section, we derive analytical expressions for the 
stationary state phonon number to ss and for the effective 
cooling rate 7 C in this equation which are much more 
general than the expressions which we obtained already 
in Section IIIIDI Afterwards, we show that these rates 
are in good agreement with numerical solutions of the 
closed set of 23 cooling equations which can be found in 
Section [niB] and App. 

A. Stationary state phonon numbers 

Using these cooling equations and setting the time 
derivatives of all expectation values equal to zero, we 
obtain the stationary state phonon number 

with the frequencies £1, £2, and £3 defined as 

g = (4A 2 + r 2 )(r 2 + z/ 2 ) + 2(r 2 + 3^)fi 2 , 
e| = (r 2 + ^ 2 ) [(r 2 + 4A 2 ) 2 + 8(r 2 -4A 2 )i/ 2 

+i6z/] + 4 [(r 2 + 2z/ 2 )(r 2 + 4A 2 ) - 81/ 4 ] n 2 

+4(r 2 + 4v 2 )n 4 , 

£f = 2(2A + /y)(r 2 + j/ 2 ) [r 2 + 4(A-^) 2 ] v 
+ [3r 4 - (4A 2 - 8Az/ - 7v 2 )T 2 
-4(A 2 - 6Av + 5v 2 )v 2 ] n 2 . (61) 

This result applies in zeroth order in rj without any ap- 
proximations and is the main result of this paper. Fig. [5] 
shows TTigg CIS Ci function of the two laser parameters f2 
and A for different experimental parameters. In the weak 
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FIG. 2: Logarithmic contour plot of the stationary state phonon number m aa in Eq. (|60[) as a function of the laser 
parameters Q and A for (a) v = 0.01 T, (b) T = v, and (c) F = 0.01 v. 



confinement regime one should choose A = 0.5 T. How- 
ever, for v ~ r and v 3> T one should choose A close 
to v in order to minimise the final kinetic energy of the 
trapped particle. 

Fig. [5] moreover shows that m ss depends only weakly 
on the Rabi frequency Q which implies that effective laser 
cooling is not restricted to laser Rabi frequencies much 
smaller than T as it is often implied in the literature 
[DL EE EE H3 ■ This is an interesting result, since larger 
O's yield higher cooling rates 7 C , as we shall see below. 
Our numerical simulations show that f2 can be as large as 
0.3 r in the weak confinement and as large as 0.3 v in the 
strong confinement regime without noticeably increasing 
the final phonon number m ss . However notice that for 
Si's larger than that, cooling changes quickly into heat- 
ing, since to ss increases relatively rapidly beyond certain 
critical points. 

Let us now have a closer look at different parameter 
regimes to emphasize that our results are consistent with 
previous calculations. For example, for relatively small 
Rabi frequencies ft, Eq. (jSH)) simplifies to 



r 4 + 8r 2 (A 2 + v 2 ) + 16(A 2 — v 2 ) 2 



(2A 



16i/A(r 2 + 4A 2 ) 
! +4(A 



4A(r 2 + 4A 2 ) 



(62) 



This equation is in good agreement with the stationary 
state phonon number implied in Refs. [HI, ED, EE HH] for 
c?3 = 0. This can be shown using Eqs. ([2]) and (|43| in this 
paper and the expressions for A± in Eq. (7) in Ref. fl9| . 
However, notice that the A± in this paper are slightly 
different from later definitions of the A± [lE E3 . 



1. Weak confinement 



In the weak confinement regime, the stationary state 
phonon number m ss in Eq. (|62[) simplifies to 



in 



weak 



16Ai 



(3r 2 - 4A 2 )» 2 
8^A 



(63) 



Exactly the same stationary state phonon number is ob- 
tained when setting the left hand side of the five effective 
cooling equations in Eq. (pLTj) equal to zero. This con- 
firms the consistency of the calculations in this paper. 
As already pointed out above, for small f2, this expres- 
sion assumes its minimum if 



A = ir. 

2 



(64) 



For this laser detuning and Rabi frequencies SI <C T, the 
stationary state phonon number simplifies to 



weak 



r 



(65) 



which is much larger than one. 



2. Strong confinement 

Using the effective cooling equation derived in Section 
IIIIDi we find that the stationary state phonon number 
in the strong confinement regimes equals 



strong 



(A-^) ; 
4i^A 3 



[(A + u) 2 - (2A + v)v\ (66) 



to a very good approximation. Exactly the same result 
is obtained when neglecting terms proportional to T and 
O 2 in Eq. (fBT)]). This result suggests immediately that 
one should choose 



A = 



(67) 



in order to minimise the final phonon number to ss . When 
substituting this detuning into Eq. (|62|) . we find that 
the stationary state phonon number for laser sideband 
cooling is for relatively small Rabi frequencies SI to a 
very good approximation given by 



16^ 2 

which is much smaller than one 



[46> - 3] 



(68) 
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FIG. 3: Logarithmic contour plot of the cooling rate 7 C in Eq. (J69J) in units of 1/r as a function of the laser parameters SI 
and A for (a) v = 0.01 Y and 77 = 0.01, (b) T = v and 77 = 0.1, and (c) Y = 0.01 v and 77 = 0.1. 



B. Effective cooling rates 



2. Strong confinement 



Proceeding as above but eliminating only the expec- 
tation values other than 712 adiabatically from the time 
evolution of the trapped particle, we find that the effec- 
tive cooling rate j c in Eq. ([55]) is in second order in 77 
given by 

16?7 2 i/ArJl 2 £f 



7c 



*2° 



(69) 



with £1 and £2 as in Eq. (|61[) . Fig. [3] shows this cooling 
rate as a function of the laser parameters ft and A. In- 
deed we find that j c increases in general as ft increases. 
This confirms that one should choose ft in general as large 
as possible without noticeably decreasing the stationary 
state phonon number m ss . 

For relatively small Rabi frequencies ft, i.e. for ft ^ V, 
the cooling rate 7 C in Eq. simplifies to 



7c = 



77 2 rn 2 



77 2 rn 2 



r 2 + 4(A 



r 2 + 4(A- 



(70) 



When comparing this result with the expression for A_ 



compai 

0E1 



we find that our cooling rate is essentially 



A, 

the same as the cooling rate reported in these earlier 
references. Moreover, as Fig. @] confirms, Eq. (|6"9")l is in 
very good agreement with numerical solutions. 



1. Weak confinement 

In the weak confinement regime, the laser detuning A 
which minimises the stationary state phonon number m ss 
is given by \ T (cf. Eq. Taking this into account, 

the cooling rate 7 C in Eq. (|69| simplifies to 



7c 



277W 

r 2 



(71) 



in the limit of weak driving, ie. when neglecting terms 
proportional to ft 2 . As Fig. @Jc) illustrates, this cooling 
rate is in good agreement with a numerical solution of 
the full set of 23 rate equations. 



In the strong confinement regime, terms which scale 
as T or ft are in general negligible (cf. Eq. (|51|0 . Taking 
this into account and simplifying Eq. (|69[) accordingly, we 
find that the cooling rate 7 C in this equation is exactly 
the same as j c in Eq. ([53")) . The cooling process becomes 
indeed the most efficient, when the detuning A is close 
to the phonon frequency v (cf. Eq. (|67"|) ). The cooling 
rate is in this case given by 



7c 



?7 2 r> 2 



(72) 



which is essentially the stationary state photon emission 
rate of a laser-driven atom multiplied by rj 2 . 



VI. CONCLUSIONS 

This paper analyses the cooling process of a single 
trapped particle with red-detuned laser light. In con- 
trast to previous authors [TJ, [H, Q, [H, HH , our analysis 
avoids the adiabatic elimination of the excited atomic 
state when calculating the effective cooling rate j c and 
the stationary state phonon number m ss . Our calcula- 
tions hence apply to a wider range of laser Rabi frequen- 
cies ft. They show that ft can be chosen relatively large 
without affecting the final outcome of the cooling process. 
As illustrated in Fig. f2J in the weak confinement regime 
with i/«r, the Rabi frequency ft can be as large as 0.3 T 
and in the strong confinement regime with rcc, the 
Rabi frequency ft can be as large as 0.3 v without reduc- 
ing the stationary state phonon number m ss noticeably. 
This is an interesting observation, since j c scales as ft 2 
and increases rapidly when ft increases. 

The main novelty of our calculations is a transfor- 
mation of the original Hamiltonian which replaces the 
atomic lowering operator er~ and the phonon annihi- 
lation operator b by two new operators x and y which 
commute with each other. The corresponding particles 
are neither atomic excitations nor phonons but provide 
a more natural description of the trapped particle. 
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FIG. 4: Logarithmic plot of the time dependence of the mean phonon number m for the experimental parameters indicated 
in the figure. Here (a) illustrates the strong confinement, (b) shows the medium confinement, and (c) shows the weak 
confinement regime. The numerical solutions are the result of a numerical integration of 23 cooling equations while the 
analytical solution represents Eq. (|59|l . 



Our calculations arc therefore more straightforward 
than previous calculations. Since the theory of laser 
cooling has already been studied in great detail in the 

literature EMI, 

the main purpose of this paper is to 
establish and test a framework for the modeling of laser 
cooling which can be extended relatively easily to more 
complex cooling scenarios like cavity-mediated laser 
cooling 0, HfjJ and possible quantum optical heating 
mechanisms in sonoluminescence experiments p5| . 
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Appendix A: Derivation of the reset state IZ(p) 

Substituting the interaction Hamiltonian Hi in 
Eq. p7|) into the expression for Tt{p) in Eq. (f2U)) and 
using first order perturbation theory we hnd that 

n(p) = lim — / dt dt'V 
At^oAtJ J 

(lkA|#lO)|0) /9(0|^i(i')|l kA ) . (Al) 

Using the concrete expression for H\{t) in Eq. (fTT)) . the 
relation 




dt' e K"*-^Xt-t') = 27rS(uJk - Uh) , (A2) 



and ignoring an overall level shift which can be absorbed 
into the free energy of the system, the reset operator IZ(p) 
becomes 

K(p) = ^27r|.g kA | 2 <5(w fe - cj l ) 

kA 

_ _ /iryfecostA ( ir]kcos'd\ , . 



This expression can be evaluated relatively easily in the 
large volume limit, where 




(A4) 

when using the same notation as in Section III Bl Per- 
forming the integration implied by this equation, using 
Eqs. (fT4|) and (fl"8|) . and writing the normalised dipolc 
moment Dqi = Doi /| Doi I as 

Doi = (d 1 ,d 2 ,d 3 f (A5) 

we finally find that IZ(p) is indeed given by Eq. ([25]). The 
spontaneous decay rate T in this equation is the same as 
T in Eq. (|24[) . Notice that d% denotes the component of 
the normalised dipolc moment Doi in the direction of the 
cooling laser. 



Appendix B: Relevant expectation values 

The calculations in the following two appendices re- 
quire in addition to the expectation values defined in 
Section MI Bl the x operator expectation values 

fci = (x + :z; t } , k 2 = i(x-x^}. (Bl) 

Moreover we employ the mixed operator expectation val- 
ues 

7i 4 = (x^xy^y) , fci3 = ((as + x^)y^y) , 
fci4 = i ((a: - xt)yiy) , k 15 = ((x - x f )(y - y 1 )} , 
ki 6 = i ({x + x r )(y - yt)} , k 17 = ({x + x r ){y + y r )) , 
kl8 =i{(x-x^(y + ^)}, (B2) 
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and 



hg 
ho 
hi 
hi 
hz 
hi 



i{( x + x i)(y 2 -yV)), 
({x + xW+yV)), 
i((x-x^(y 2 + y^)}, 
(xix(y 2 + yV)), 
i(Jx(y 2 -yV)). 



(B3) 



The time derivatives of these and other expectation val- 
ues which we defined in Sec. lIII Bl can be found in App.lCl 
Notice that this paper does not introduce expectation 
values 77.3 and £3 to h. These names are reserved for 
variables which are used in other cooling scenarios. 



Appendix C: m, fcn, and ki2 in the weak 
confinement regime 



Setting 77 = and substituting the x operator expec- 
tation values ni, k%, and h into Eq. (|36[) . we find that 
they evolve according to 



1 



"1 

h 
h 



-Q h — r i%i j 



1. 



-rfci 



- 2ni) + Afci - -Tk 2 



(CI) 



in zeroth order in 77. These equations form a closed set of 
differential equations. Eliminating all x-operator expec- 
tation values adiabatically which change on the relatively 
fast time scale given by T and adopting a notation where 
x = .t' ' + x^ + with the superscript indicating 

the scaling of the respective term with respect to 77, we 
find for example that n\ is in zeroth order in 77 given by 



and hi and ki5 to hs is in zeroth order in 77 given by 
1 



fell 
ki2 
ki5 
kw 
kn 
km 



-Clk 



18 



vk\2 — rfcn 



1 



:^hi 



-Qfci5 + vkn - Tki2 , 
-fi(fc 8 - 2/C12) - A he — v k ls 

Afci 5 +L>ktf- ^ rfc 16 > 
-Afcig -1//C16 - ^ rfc 17 ! 



Q(k 7 - 2ku) + A kir + v fciB ~ 2 F fcl8 ' ( C4 ) 



All six expectation values evolve on the relatively fast 
time scale given by the spontaneous decay rate T. Taking 
this into account and eliminating them adiabatically in 
the weak coupling regime, i.e. for relatively small v, we 
find that 



k i°i = ^ i^h - (3r 2 - 4A>fc 8 ] , 



n 



-(0) 
' 12 



M 4 r 
n 2 
M 4 r 



[(3r 2 - 4A V fc 7 + /j 2 r fc 8 ] (C5) 



to a very good approximation. The constant /z 2 is given 
in Eq. (|C3j) above. In addition to and k{^ we obtain 



expressions for fcjg' and fejg . These arc used in the next 
subsection to calculate n% up to first order in 77. 

Proceeding as above but taking terms up to first order 
in 77 into account we find that the first order in 77 contri- 
butions of the x operator expectation values m, h, and 
h in Eq. (|4Q]) evolve according to 



'.(1) 



1 



nk 



(i) 



r n 



(1) 



-r\v k 



.(0) 



Ak 



(i) 



-Tk 



(i) 



k « = -2f2n« - ^ fc$ + A jfeW - -T k^ 



(C6) 



,(°) 



ft* 



The constant /z 2 in this equation is given by 



f i z = 2n 2 + r 2 + 4A 



(C2) 



(C3) 



These equations form a closed set of cooling equations, 
when the results for kf^ and k^ which we obtained in 
App. [Cl are taken into account. Eliminating m, h and 
h adiabatically 



8??z/Aft 2 



(C7) 



JO) 



In addition to n^°\ we obtain solutions for k[ and k 
These will be used in the next subsection to calculate the 
coherences hi and ki2 up to first order in 77. 

Setting 77 = and using again Eq. (|36|) . we find that 
the time evolution of the mixed operator coherences hi 



in the weak confinement regime which we introduced in 
Section MI Cl This means terms proportional to t^ 2 have 
been neglected. 

In order to calculate hi and fci2 up to first order in 
77, we need a closed set of cooling equations which holds 
correctly up to this order. Applying Eq. (|36[) again to 
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fen and ki2 and ki5 to k%g, we find that 



k" 



1 



ni 2 

(i) _ 

l 12 — 



.(1) 



,(0) 



(1) 



-Q, k[^ + v k\± — T k[^ , 



K 15 



2fl k[l ] - A k^' - v kY 8 > + rjv { k\ l " + 2k\" 3 



(i) 



(i) 



,(o) 



,(0) 



"■21 ) 2 15 ' 



(o) j_ oi.(o) _ fc (o) N 



fc^ — Afe^ ^ ^16 ; + V 1 ' J 2^^i 7 

K 18 



irfc (1) 
2 16 ' 



JO) ,(0) 



(1) 



-20 Jfc£> + A fc« + v fc^ + w [kf - jfeg?) 



"2 rfe i8 ) - 



(C8) 



Substituting the definitions of the mixed-particle expec- 
tation values 714, fci3 and kn and kig to &22 into Eq. (|36[) 
and setting 77 = 0, we find that 



while 



fi4 = —£l fci4 — r ri4 , 



fel 3 



-Afci4 - 2 rfc i3 



fei4 = fi(n 2 -2n 4 ) + AA:i3--rfci4, (C9) 




FIG. 5: A comparison of the analytical results for ni, ku, and 
ki2 in Eqs. ijct]) . (fC5|) . (fC7| . and (|C11|I with the results of a 
numerical solution of the above cooling equations for r\ = 0.1, 
O = ^ = 0.01 T, and A = 0.5 Y. 



result of a numerical solution of the above cooling equa- 
tions. Very good agreement between both solutions is 
found which suggests that the effective cooling equations 
in Eq. (|45[) apply after a very short transition time of the 
order 1/r. 



fci 9 = -O(fcio - 2fo4) - A fc 2 o - 2v k 2 2 - k 19 , 
fc 2 o = Afcig + 2zyfc 2 i - 2" rfc 20 ' 



hi 



-Afc 22 - 2^fc 20 - 2 rfc21 



fc 2 2 = Q(k 9 - 2fc 2 3) + A fc 2 i + 2v fcig - -T fc 22 
^23 = k 22 - 2v k 24 - T k 23 , 
fc 2 4 = —Slkig + 2vk 2 ?, - r/c 2 4 • 



(CIO) 



These final six differential equations hold in zcroth order 
in 77. Setting the right hand side of these and of the 
cooling equations in Eq. (|C8|) equal to zero, we finally 
find that 

*S } = [2A* 10 + r], 



(i) _ 8771/Afi 2 

c 12 



[2n 2 - fc 9 + 1] (Cll) 



in the weak confinement regime. This means, terms of 
order v 2 have again been neglected. Fig. [5] compares 
the above analytical results for n\, ku, and ki 2 with the 



Appendix D: m, fcn and ki2 in the strong 
confinement regime 

The calculation of n\ in zeroth order in 77 is the same as 
in App. O However, in the strong confinement regime, 
the expression in Eq. (|C2[) simplifies to 



,(0) 



4A 2 



(Dl) 



Setting 77 ~ and eliminating the ?/ operator coherences 
adiabatically from the system dynamics we immediately 
find that k-j to fcio all equal zero in zeroth order in 77, 



„(°) 



.(0) 



JO) 



„(°) 



0. 



(D2) 



Taking this into account when eliminating the mixed 
operator expectation values whose time derivatives are 
given in Eq. (|C4[) . we moreover find that 



.(0) 



,(0) 



0. 



(D3) 



To calculate the coherences fen and fci2 up to first order 
in 77, we have a look at the time derivatives of fen, fei2, 
and fei5 to fcis in first order in 77 which can be found in 
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FIG. 6: A comparison of the analytical results for m, fen, 
and feia in Eqs. (|D1|> . (|D3|) . and (|D5|l with the results of a 
numerical solution of the above cooling equations for 77 = 0.01, 
n = r = 0.01 and A = v. 



Eq. (|C8l) . Combining Eqs. (|C10|) and (|D2|) . we immedi- 
ately see that 



do) 

^19 



.(0) 



JO) 



.(0) 



0. 



(0) JO) ,(0) 
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Taking this and the expressions for k\ 
k\2 obtained in App. [C] into account, when setting the 
time derivatives of the relatively fast evolving variables 
in Eq. (|C8|) equal to zero, we therefore find that 
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These coherences arc different from the coherences in 
Eq. (|C11|) . since they apply only in the strong confine- 
ment regime. As Fig. [5] shows there is again very good 
agreement between the analytical and the numerical so- 
lutions for ni, fell; and ki2- This means that the effec- 
tive cooling equation for the strong confinement regime 
in Eq. (|52")) too applies after a relatively short transition 
time. 
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